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Abstract 

The main formalisms of partial level densities (PLD) used in preequilibrium nuclear reac- 
tion models, based on the equidistant spacing model (ESM), are considered. A collection of 
FORTRAN?? functions for PLD calculation by using 14 formalisms for the related partial-state 
densities is provided and 28 sample cases (?3 versions) are described. The results are given in 
graphic form too. Composite (recommended) formulas, which include the optional use of various 
corrections, i.e. the advanced pairing and shell correction in addition to the Pauli effect, and 
average energy-dependent single-particle level densities for the excited particles and holes, are also 
given. The formalism comprises the density of particle-hole bound states, and the effects of an 
exact correction for the Pauli-exclusion principle are considered. 
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PROGRAM SUMMARY 

Title of program: PLD 
Catalogue identifier: 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland 

Licensing provisions: none 

Computer for which the program is designed and others on which it is operable: PCs (486 and 
Pentium) 

Operating systems under which the program has been tested: DOS 
Programming language used: FORTRAN 77 (MS-FORTRAN v.5.0) 
Memory required to execute with typical data: 491 Kbytes 
No. of bytes in distributed program, including test data, etc.: 1,256,659 

Distribution format: ASCII 

Keywords: Partial nuclear level density, nuclear level density, single-particle level density, 
equidistant-spacing model, preequilibrium emission, nuclear reactions 

Nature of physical problem: 

This Fortran code is a collection of subroutines for calculation of the partial nuclear level densities 
(PLD) mainly used in preequilibrium nuclear reaction models, by using 14 formalisms for the 
related partial state densities (PSD). 

Method of solution: 

The main approaches to the calculation of the partial state density, based on the equidistant 
spacing model (ESM), are used. Composite (recommended) formulas including optionally various 
corrections, i.e. the advanced pairing and shell correction in addition to the Pauli effect, and 
average energy-dependent single-particle level (s.p.l.) densities for the excited particles and holes, 
are also involved. The density of the particle-hole bound states is moreover comprised, and the 
effects of an exact correction for the Pauli-exclusion principle are considered. 

Restrictions on the complexity of the problem: 

Although the quantum-mechanical s.p.l. density and the continuum effect can also be reproduced 

by a corresponding Fermi-gas formula, to be used accordingly within the average energy-dependent 
PSD in multistep reaction models, this effect is not included. The calculation of PLD with linear 
momentum, of first interest for modelling preequilibrium-emission angular distributions, is not 
available either. 

Typical running time: 

The execution time is strongly problem-dependent: it is roughly proportional to both the number 
of the excitation energies and the exciton configurations considered in the calculation, while con- 
sistent differences arise when various PSD formalisms are used. Twenty-eight sample cases with 
73 versions which require from 0.1 to 1661 s on a PC Pentium/166MHz are provided. 

Unusual features of the program: 

The PSD functions have been optimized for their independent use, in order to provide tools for 
PSD/PLD users (see [5]). The related drawback is the increase in execution time, while a proper 
use would involve the calculation of some coefficients only once in the main program. Second, the 
PLD. FOR has been organized so that various formulas and versions may be tried as well as the 
comparison between their predictions. 
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LONG WRITE-UP 



1. Introduction 



The particle-hole excitations caused by the nuclear reactions which proceed through a num- 
ber of nucleon-nucleon interactions are described within either the semiclassical models or the 
quantum-statistical theories of the preequilibrium emission (PE) [|l],§ by means of the particle- 
hole state densities. Basic approaches to the partial state density (PSD) consist in combinatorial 
calculations performed in the space of realistic shell-model single-particle levels (s.p.l.) [^. Lenske 
et al. 1^ have already used them in order to connect in a consistent way the quantum-statistical 
theories of the multistep-direct (MSB) and multistep-compound (MSG) processes 0. However, 
the strong dependence on the s.p.l. basic set is the main of several shortcomings inherent in the 
method (e.g., 0-^] and references therein). 

The equidistant spacing model (ESM) state density including the effect of the Pauli- 
exclusion principle |11| is still widely used, as well as the phenomenological s.p.l. density value 
g ~ ^/14 MeV~^. Basic developments of the Williams formula fll]] are due to Betak and Dobes 



p!2| including the nuclear-potential finite depth correction, Stankiewicz et al. [|T3| and Oblozinsky 
TM who added the bound-state condition, and Fu HTSll and Kalbach [ITBl who included an ad- 



vanced pairing correction. Additional studies along this line have involved exact Pauli-correction 
calculation [p"^^ -[T9[|. Kalbach pO|-^ also discussed different energy dependences of the s.p.l. spac- 
ings and pointed out the necessity to study this subject closely related to PE surface effects, 
due to the interdependence of the corresponding assumptions. PSD including different energy- 
dependences of the excited-particle and hole state densities has recently been used |^ in the 
geometry-dependent hybrid (GDH) model [0,^. A similar attempt [^,0 has focused on the 
MSD and MSG processes in the framework of the Feshbach-Kerman-Koonin (FKK) theory 
An independent semiclassical analysis I^S],^ has additionally justified the surface localization of 
the most important first nucleon-nucleon interaction within PE processes, and provided average 



quantities useful for the corresponding PSD calculation ||3^ . 

The various PLD formalisms based on the ESM Williams-type formula, which are still ex- 
tensively used in nuclear-reaction calculations, determined the need for an appropriate subroutine 
collection. It is on the request of a project concerning a reference input-parameter library for 
nuclear model calculations |^ that the present work is based. Thus, the program PLD. FOR is a 
collection of algorithms developed until now and widely used for PSD/PLD calculations. The one- 
and two-fermion system versions of six different approaches and one composite (recommended) 
formula including various corrections are available as FORTRAN77 functions. 

The main points of the PSD and PLD formalisms are presented in Section 2. At the same 
time, the sample cases (Table I) for the program PLD. FOR are described. The structure of the 
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program and the input-data description is given in Section 3. Finally, a worked example is pre- 
sented in Section 4. 



2. Formalism 

2.1. Partial state density in the uniform spacing model 

The state density of a system of p excited particles above the Fermi level and h holes below 
it, considered within the uniform spacing model (based on a constant spacing d=l/g between the 
non-degenerate single-particle levels) at the last occupied level in the ground state of the nucleus, 
was obtained by Williams 



.(P./».B) = ^^||-^. (1) 

p\h\[n — Ij! 

by decreasing the excitation energy E with the correction for the Pauli blocking 

^ ^ p{p+l) + h{h-l) _ h_ 

Ag 2g ^ ^ 



with respect to Ericson early formula for n=p + h total number of excitons. The sum of the 



partial state densities for all allowed particle-hole numbers p=h is consistent with the total nuclear 
state density formula obtained in the frame of the ESM of the one-component Fermi gas 



exp[2(fgi;)'/^] 



The asymptotic equality 



uj^{E)c^Y.^iP^h,E), (4) 



is illustrated in Fig. 1(a) (the sample cases 1/lA) as proved by Williams for the generic value g=l 
MeV~^ of the s.p.l. density. 

Similarly, in the case of the two kinds of fermions considered, with the g-,^ and g^ being 
the single-proton and single-neutron state densities, respectively, the PSD for and h.,, proton 
particle and hole numbers, respectively, and py and neutron particle and hole numbers {n = 
P-K + K + Pu + K) is 

Uip.,p., K, K, E) = , (5) 



where the Pauli effect correction has now the form 



-3 



p-k{p-k + 1) + K{K - 1) ^ Pujpu + 1) + KjK - 1) 

97T Qu 



2 VS'tt Qv 



The corresponding ESM total nuclear state density for a two-fermion system with an average 
total single-particle state density g=gn+gu and related level density parameter a={7i^/6)g 

y^ exp[2(aE)V^] 
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is also consistent with the sum of the partial state densities (^) for all allowed pairs of particle-hole 
numbers Pn=hT, and p^=h^, i.e. 

u;2(^) ~ uj{Pn,Pu,K,K,E) , (8) 

shown in Fig. 1(b) (Cases 2/2A) for the similar generic values g.^=gy=g/2=l MeV~^. It results 
that Eq. (P) is true within 4% for E > 3 MeV. A comparison of the PSD given by the Williams 
one- and two-fermion formulas for the real case of the nucleus ^^Nb is carried out in Fig. 1(c) 
(Cases 3/3A) by using the phenomenological value g=A/13 MeV~^ and the derived quantities 

9^ = -j9 (9a) 

9u = -^9 ■ (9b) 

The renormalization of the PSDs given by one-component Williams-type formula can be done by 
using the ratio between the total state densities given by the two- and respectively one-component 
Fermi gas formulas of the general form [0^ , ^ 



_ exp|2K0!5 (10) 



v/? exp[2(a^)V^] 
'^^^^^ = T2"aV4(f/ + t)5/4 ' (11) 

where the nuclear temperature t corresponding to the effective excitation energy f/ (see below) is 
defined by 

U = at^ -t. (12) 
Thus, the renormalized one-fermion PSD has been defined by 

u^{n,E) = j{U)uj^{p,h,E), (13) 
where the renormalization factor j{U) (two-fermion system correction - TFC) 



7r 



1/2 U 



has only a weak energy dependence approximately equal to [/~^/^. 

Since each of the closed formulas (|iy) and ([TT|) are asymptotically equal to the sum of the 
one- and respectively two-component PSDs over all allowed particle-hole numbers p=/i, the corre- 
sponding sum of the renormalized PSD (0) is consistent with the closed formula for a two-fermion 
system, as shown in Fig. 1(d) (Cases 3/3A/3B). 

2.2. Bound-state and finite depth corrections 

The limitation on the hole maximum energy due to the finite depth of the nuclear potential 
12| as well as on the maximum particle excitation by the nucleon binding energy B in the case of 



the bound states, yielded the approximate one-fermion ESM formula [in 
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u{p, h,E) 



p\h\{n 



(-1)^^' C; CUE - A,, -iB- jFr-'e{E - a,, - iB ~ jF) 



(15) 



where F is the Fermi energy which is now considered to be halfway between the last filled and the 
first free s.p.l. in the nucleus ground state in order to have a PSD form symmetrical in p 

and h. Under these circumstances, the Pauli correction term in Eq. ( [T3| ) is 



A 



ph 



p{p - 1) + h{h - 1) 
'^9 



while 



ph 



P^ + /l^ 

2^? 



(16) 



(17) 



is the minimum energy of the (p, h) state due to the Pauli blocking. 6 in Eq. ( |T5| ) is the unit step 
function, i.e. 1 for a positive argument and otherwise. The effect of the nuclear-potential finite 
depth on the PSD as firstly pointed out by Betak and Dobes |T2|, and the additional one due to 
the condition of bound states, included by Oblozinsky [|14|, are shown in Fig. 2 (Case 4) and Fig. 
3(a) (Cases 5/5A/5B). The changes obtained by releasing consecutively the finite-depth potential 
and bound state conditions, by means of large F- and -B-values, are shown in the later case for 
the basic exciton configuration Iplh. The consequences of these conditions on the total nuclear 
state density, as well as the corresponding results obtained by using the asymptotic formula (H), 
are illustrated in Fig. 3(b) (Cases 6/6A/6B). 

The same method [|1^] applied within the two-fermion ESM led to the density of partial 
two-fermion bound states of the form 



u{p^,K,p^,h^,E) 



pPw + hn pPv + hiy 



Ptv Pi/ Ht, hi, 

E\ " \^ \^ /I YTT+iv+iiT+jv ri^TT r-'ii' r^^^ /^•J" 
U 'IL 'U 'li 'Vli- LV ^ ^ ^> ^P^^P.^h^^h, 



^ {^E Ap^fi^p^fi^ ij^Bj^ iyBy jirF-jj jpFv 



\n-l 



where 



PtV i^TT -iPu ihi/ 



pAptt - 1) + K{K - 1) _^ Pujpu - 1) + K{K - 1) 



9^ 



9u 



and 



_ pI + K ^pI + K 

Oip^^h-^^p^^hi, ^ + 



29. 



29u 



(19) 



(20) 



are symmetrical in the respective particle and hole numbers. The total nuclear state density, i.e. 
the sum of all PSD given by Eq. (^ for allowed pairs of proton and neutron particle- hole numbers, 
calculated with/without the potential finite-depth and bound-state conditions are shown in Fig. 
3(c) (Cases 7/7A/7B). 
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2. 3. The advanced pairing correction 



2.3.1. One- component Fermi gas formula 

(a). In order to take into account the pairing interaction between nucleons a correction was 
included [|15| into the formula derived by Williams in the frame of the free Fermi-gas model (FGM), 
function of the particle and hole numbers as well as excitation energy of the configuration. The 
Pauli correction was also modified to be consistent with the pairing correction, so that the PSD 
formula became 

u^{p, h,E,P + B) = , (21) 

p\h\[n — 1)! 

where B is the modified Pauli correction following the Williams term A: 

B = A[l + {2g^/nfY'^ , (22) 

and the pairing correction term 

P = \g{Al - A') (23) 

is determined by the ground- and excited-state gaps Aq and A{p,h,E). The former is related 
to the condensation energy C = gA^/A. On the other hand, the constant pairing correction Up 
for the total state density, based on the odd-even mass differences (e.g., [^), may be rather well 
related to the value P(n) |T3| where n is the most probable exciton number. Since A=0 if n > n, 
it results that Aq can be derived from the relation Up = gAl/A. Then, Fu obtained the following 
parametrizations for A ||15|JT6[| : 



A = 0.996 - 1.76(nK)i-'^°(E/C)-°-''^ if E > Ephase 

^0 

= if E< Ephase , (24) 

where ric = 0. 792(7 Aq is the critical number of excitons and Ephase is the energy of the pairing 
phase transition given by 

Ephase = C [0.716 + 2.44(n/n,)2-i7] jf ^/^^ > o.446 

= if n/n^ < 0.446 . (25) 



Actually, the lower limit in Eq. (^5]) was adopted |]T6l in order to take into account explicitly the 



lack of a phase transition for small n. The original comparison of the calculated PSD with and 
without pairing correction is shown in Fig. 4(a) (Case 8, with the threshold-condition released). 
One should note that the respective equations provide PSD- values even below the minimum exci- 
tation energies (thresholds) characteristic of each configuration 

Uth = C [3.23(ri/ne) - 1.57{n/n^Y] if n/n^ < 0.446 

= C [1 + 0.627{n/ncf] if n/n^ > 0.446 . (26) 

The correct densities, i.e. above the respective thresholds, are shown in Fig. 4(b) (Case 8). 

The total state density obtained as the sum of all PSD for allowed pairs of particle-hole 
numbers p = h are compared with the closed formula of the one-component Fermi gas used at an 



effective excitation energy decreased by the constant pairing correction Up |T^: 
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(27) 



A distinct underestimation by this formula of the sum of PSD for E <8-10 MeV also results from 
the comparison shown in Fig. 4(b) (Case 8) for the constant pairing correction Up=3.5 MeV. The 
analysis performed with respect to the effect of various Up values on this agreement |[T^ is shown 
in Fig. 4(c) (Case 9) for the generic values g=4: MeV~^ and 0, 2, and 4 MeV for Up. 

The relationship between the pairing corrections for the PSD and the ESM total state density 
was extended to considering the nuclear-shell effects by using an additional back-shift S of the 
effective excitation energy |15| 

u{p, h,E,p+B+s)= ^"^^ - - . 

p\h\{n — 1)! 

By analogy with the BSFG model for the two-fermion system total state density 

y?exp{2[a(E- A)]i/2} 



(2^ 



12 aV4 (E - A + t)5/4 



(29) 



the back-shift parameter is connected to the BSFG virtual ground-state shift parameter A through 
the relation 



A 



Up + S 



(30) 



Thus, the predictions of the closed formula of the one-component Fermi gas (|27D for the effective 
excitation energy U = E — A, and the total state density standing for the PSDs given by Eq. 
(|28|) , should be asymptotically equal as shown by the upper curves in Fig. 4(d) (Cases 10/lOA). 
The lack of consistency between the sum of the one-component PSDs and the two-fermion BSFG 
state density formula has been illustrated by means of the various predictions of one- and two- 
component closed formulas with the same parameters for the excited odd- A nucleus ^^Ca ||T^, also 
shown in Fig. 4(d). 

(b). An improved implementation of the advanced pairing correction within the Williams 
formula 



u{p,h,E,AK) 



g^[E-AK{p,h)] 



(31) 



p\h\{n - 1)! 

adopted a Pauli correction function symmetric in particles and holes, which also included the 



effects of a passive hole |TB 



Ak{p, h) = Ethreship, h) 



p{p+l) + h{h + l) 
^9 



(32) 



where Ethresh{p, h) = pf^/g and pm = maximum{p, h). The inclusion of the pairing interaction 
led to the modified form of the threshold energy for a given exciton configuration 



E, 



thresh 



{p,h) 



4 



Pn 



Prn\ ^^2 



1/2 



(33) 



A third term was added to the modified Pauli-and-pairing correction 



Ak{p, h) = Ethreship, h) 



p{p + I) + h{h + I) , {p-lf + {h-l) 



45 



+ 



gF{p, h) 



(34) 
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where 



F{p, h) = 12 + 4g[E - Ethreship, h)]/p„ 



(35) 



in order to obtain the PSD values of g and 2g for E = Ethresh{p, h) and E = Ethresh{p, h) + 1/g, 
respectively. The comparison of the values given by the formulas of Fu []T3| and Kalbach |TB[ above 
Ethreship, h) for cach exciton configuration, respectively, is reproduced in Fig. 5(a) (Case 12). 

The pairing effects were also included in the total state-density formula by using an effective 
excitation energy decreased by the effective pairing shift |jl6[ 

PeffiE, C) = maximum{E2, C/{1 + exp[4(0.625 - E/C)]}) , (36) 

where the quantity E2 is written as 



E2 = C[1 + 2.508/(nj2] 
= C [6A6/nc - 6.28/(nc)2] 



if fic < 4.4S 
if rir > 



(37) 



The sum of the PSD provided by Eq. (^) is now consistent with the Fermi-gas formula (|37|) 
provided that the constant Up is replaced by the effective-energy shift Peff{E, C), as shown in Fig. 
5(b) (Case 13). 

(c). The PSD formula also within the ESM and based on an exact calculation of the Pauli 
correction term [|TB| , extended to the case of the finite well depth and bound states, and including 
the Kalbach |jl6[ pairing correction, is |[T9[] 

P h 



uj{p,h,E) 
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n-l 



^,uT.T.i-^y^'c;ciY.{E-E-, 

P-"'- i=o j=0 A=0 



thresh 



iB-jF) 



n-l-A 



X e(E - Ethresh -iB- jF)B{p, h, A)- , (38) 

[n — 1 — Ay. 



where the coefficients B{p, h, A) have the expression 



B{p, h,\)=Y: C{p, Xi)C{h, A - AO , (39) 

Ai 



and the coefficients C(m, A) are determined by the recursive relation 

C(m, A) = ^k{-m/gyC{m - 1, A - z) , (40) 

with 



C(0, A) = 1 for A = 

= for A 7^ , (41) 

while the Bernouli numbers bi can be found in mathematical tables. Eq. (|38| ) becomes the PSD 
expression of Baguer [l^ in the limiting case of large B and F. The exact coefficient B[p, h, A) 
corresponds to the factor 
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which follows from the approximate formula of Oblozinsky [|T^] if the energy-term power is ex- 
panded by means of the binomial theorem. The comparison with the PSD values obtained by 
means of Kalbach formula and parameters g=14: MeV~^ and Ao=l MeV, shown in Fig. 5(c) 
(Case 17A), proves the suitability of the Pauli-correction approximation. At the same time, the 
total state density given by the sum of the corresponding PSDs calculated by using the parameter 
global values g=8 MeV~^ and Ao=l MeV [p^p^ is in good agreement with the one-fermion closed 
formula (|42|) by using the effective-energy shift Peff{E,C) [Fig. 5(d) and Case 17]. 

The comparison of Eq. (|38|) with the Oblozinsky formula can be carried out by neglecting 
the pairing effect, i.e. through substitution of the threshold energy Efhresh by the Pauli energy 
aph- While Oblozinsky results for small exciton numbers are well represented (Case 15), a relative 
deviation of about 10-40% exists between the two sets of calculated PSDs when larger 2p3h and 
3p2h exciton numbers are involved [Fig. 6(a) and Case 15A]. On the other hand, the inclusion of 
the pairing effect decreases the partial bound-state densities for small E and enhances them 
for large E [Fig. 6(b) and Case 16]. 

2.3.2. Two- component Fermi gas formula 

(a). The PSD formula derived by Williams in the frame of the two-component free Fermi-gas 
model, became by inclusion of the pairing correction P2 |p6| 



o;(p., K,p.,K, E, P, + B, + S)= (I) " ,, , (43) 

\2J pjhjpjhj{n-iy. 

where B2 is the simple extension of the one-fermion system correction factor for the Pauli-exclusion 
principle modified by the pairing effects. The pairing correction term has been adopted under the 
assumption of no pairing interaction between the protons and neutrons so that 

P2(^,n^,n^) = PiiE^,n^) + PiiE^,n^) , (44) 

where 

Pi{E^, n^) = ^g^Al^ - Al{E^, n^)] , (45a) 

Pi{E,,n,) = ^g,[Al - Al{E,,n,)] , (45b) 

and E = Et, + Ey, tIt, = Pt, + Ht,, and n^, = p^ + h^. Based on the pairing theory for two kinds 
of fermions, an approximate solution was adopted for the gaps A^^ and A^. Actually, by using 
the mean gap-approximation, i.e. gT,=gy and Ao7r=Ao,y, it was shown that the results obtained for 
the one-fermion system can be used for each of the two systems, while the pairing theory yields 
that the proton system and the neutron system are excited isothermally [^. Next, the following 



simple procedure was chosen to define E^^ and E^, to approximately 10% of the exact values except 
for energies near the threshold 

Et, = Enj^/n , (46a) 
Ey = Enjn , (46b) 

with the overall error in P2 estimated to be about 2% except for energies near the threshold. 
Moreover, similar to the one-fermion case, the pairing correction Up for the total state density is 



related to P2(-£', "/^TT, "/^i/) ||36 



Up = P2{E,n^,n,) , (47) 



10 



so that 



A^, = = AUjg , (48) 

i.e. the ground-state pairing gaps in the proton system and the neutron system are both equal to 
that of the one-fermion system if the same values of g and Up are used. 

(b). The improved implementation of the pairing correction given by Kalbach |TB| has the 
two- component version 

. > u w A ^ 9n^9u'[E - AKjPn, K,Pu, KT'^ ...^ 
.{p., K, p., E, A,) = pjhjpjhjin-l)l ' 

where 

Ak{P7,, K,Pu, K) = AxiPn, K) + Ak{Pu, K) , (50) 

and the respective one-fermion expressions are used except the function F{p, h) which now has 
the form [ITB 



F{pi, hi) = Urii/n + Agi[Ei - Ethresh{Pu hi)]/prm , (51) 

where i can be either vr or u. The average excitation energies of the two kinds of nucleons have 
also been assumed to be proportional to the number of degrees of freedom of each type. 

Additional comments should concern Kalbach 0] energy- dependent pairing corrections for 



the consistency between the total state-density closed formula of the two-component Fermi gas 
at the effective excitation energy decreased by the effective-energy shift Peff{E,C), and the sum 
of all PSD provided by Eq. (|^) for allowed pairs of particle-hole numbers p,r=^7r and Pi,=h^, 
as shown in Fig. 7(a) (Case 13A). It was expected that the two-fermion PSDs sum should be 
well approximated by the Fermi gas formula when P=C.,r+Ciy, at excitation energies not too close 
to threshold, under the assumption of the same gap parameters Aott and Aqu- However, the use 
of the sum Peff^E^C-,,) and Peff{E,Cy) as well as of the constant P=Up is shown in Fig. 7(b) 
(Case 13A) to be correct only for E/C >2. The PSDs sum is underestimated, i.e. the pairing 
correction is overestimated at the lowest energies where the exciton configurations of only one 
kind of fermions are significant [Fig. 7(a)]. This suggests the use at these energies of only one 
correction term out of the two Pe//(-E', Cvr) and Peff{E,Cy). The gradual inclusion of the second 
one, by means of the sum Peff{E,C.,r)+xPeff{E,C,y) with x between and 1 for E/C varying 
from 1 to 2 as also shown in Fig. 7(c), is correct just above the first threshold. Therefore, we 
found the following form was necessary to obtain the closed-formula predictions rather consistent 
with the PSDs sum as in Fig. 7(b) 

P^+.,eff{E, C) = P^ffiE, a) (^1^ + l) , (52) 

where i can be either tt, or z/ or an average of the two terms, while the quantity E2 

E2 = Ci[l + 2.508/{nc^f] if < 4.48 

= Q [6.46/n,, - 6.28/(n,,)2] if n,, > 4.48 , (53) 



is the ground-state threshold energy for the nj=2 states, which Kalbach [|T^ also involved in the 
definition of Peff{E,C). The reduced consistency of the two curves just above E=C could be 
shifted to the region above the threshold by replacing the energy E2 with the zero value, i.e. by 
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carrying out the transition described by Eq. ( [5^ ) in the energy range between and C. The 
corresponding pairing correction is denoted by P^^^^^j in Fig. 7(b). 



The comparison of the two-fermion PSDs sum with the sum of the one-fermion PSDs [ITB] 



multiphed by the two-fermion system correction |32,34| for all allowed particle- hole numbers p = h 



as well as of the corresponding closed- formula values, is shown in Fig. 7(d) (Cases 13A/13B). The 
agreement of the related quantities is quite good but only at E/C >2. The use of Peff{E,C) 
within the two-component Fermi-gas total state density, which is the case when the two-fermion 
system correction is used, leads once again to an underestimation around the condensation energy 
C. This could be the main limit of the TFC method, which reveals the need for a different pairing 
correction at energies E/C <2 within the complete two-fermion system approach. 

(c). The PSD formula within the ESM and exact calculation of the Pauli correction term 
13 , extended to the case of the finite well depth and bound states and including the Kalbach [|TH| 



pairing correction, has the following form in the two-fermion system case 

„nn qTli, Pn fe-TT hi, 

u{p.,K,p.,K,E)= E E E E(-ir^^'"^^^^"^;:^r.c^;:c^r. 

PTT-lt-TT-Pu-'t-U-yt- J-J- i^=0 jn=0 1^=0 j^=0 

n— 1 1 

E r-'-'9it)A{p^,K,p,,K,\)- — , (54) 

t^o (n-1 -A)! 



where 



t — E Efj.f,sfi ^TT-S-TT jnFn iyBi, jyFy , (55) 



A{p^, K,Pu, K, A) = E ^(^'f' K)B{p^, X- \^) , (56) 

A^=0 



while the coefficients B{p,h,X) are determined by the Eqs.( ^9| - |41|) . The total state density ob- 
tained as the sum of the PSD for all allowed particle-hole numbers is shown in Fig. 8(a) (Case 
17B), as well as the comparison with Kalbach closed formula including the above-discussed ef- 
fective pairing correction, the only constant pairing-correction Up |l^, and the energy- dependent 
correction P^//, shown in Figs. 8(b) (Cases 17B/17C) and 8(c) (Cases 17A/17D). Fig. 8(d) (Cases 
17A/17D) also gives the comparison with the corresponding PSD values obtained by means of the 
Kalbach formula |]TB[ , which supports the correctness of the approximation for the Pauli-correction 



term. 

2.4- Partial state density with surface effects and energy- dependent s.p.l. densities 

The surface effects which may be considered within the initial target-projectile interaction 
22| by means of the PSD formula introduced by Kalbach [^, 



uip, h, E, F) = u;{p, h, E, oo)/(p, h, E, F) , (57) 



where the finite-depth correction in addition to Lu{p, h, E, oo) given by, e.g., Eq. (pTj), is brought 
off by the function 



fiP, h, E, F) = E (-1)^' Ci (^^P^^ e(E - jF{h)) 



{hi 
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where the Fermi energy F=38 MeV corresponds to the central nuclear well. The shallower potential 
in the region of the nuclear surface, where the first target-projectile interaction in PE reactions is 
most probably localized (see also p7| , p9| ), is taken into account by using an average effective well 



depth and the related Fermi energy Fi for the hole number h < 2 

F{h) =Fi for /i < 2 (59a) 

= Fo = 38MeV for /i > 2 , (59b) 

so that Fi-value is taken into account for only the initial configurations, following the assumption 
of the surface localization of PE two-body interactions exciting them. 

The smaller effective well depth may also determine an increased significance of another 
effect, namely the s.p.l. energy-dependence of g{e). First, by taking into account this dependence, 
separate excited-particle state density gp and single-hole state density gh are involved. Second, the 
increase in gp has generally compensated for the decrease in gt except when the reduced potential- 
well depth makes the latter much closer to the value at the Fermi energy, gQ = g{F); since the 
excited particles may well be excited above the Fermi level in PE reactions, opposite to the case 
of the statistical equilibrium, the gp could increase significantly. 

The interdependence of the PE surface effects and the energy dependence of the s.p.l. density 
makes the functional form of the latter to be yet an open question p2|. Thus, the first-order 



effects of the FGM dependence have been adopted P^J33J23] for the time being 



I £V'^?4f£)- (60) 



FJ 2F \FJ 

Next, the actual basic point consists in the use of the average excitation energies Up = Ep — F for 
excited particles, and Uh = F — Sh for holes [^. The former has been estimated to first order 
from the ESM Eq. (pTl), with the result 

_ E f{p + l,h,E,F) 

"^"n f{p,h,E,F) ' ^^'^ 

while the related one for holes is 

E -pUp 

Uh = , (62) 

both of them reducing to E/n if the finite depth of the potential well is not considered. The 
corresponding average FGM s.p.l. densities for the excited particles and holes 



gp{up) = go — = gp{p, h) (63a) 



Fo 


+ Up 




Fo 


Fo 


- Uh 


Fo 



1/2 

gh{uh) = go i^^^^jr^j = gh{p, h) (63b) 
are used within the ESM partial state-density ( |57D which becomes 

p\h\{n — 1)\ 

It should be underlined that the effective well depth Fi is involved in evaluating only the quantities 
Up, Uh and /(p, h, E, F) while the central well depth Fq is assumed in calculating the effective 



s.p.l. densities. The factor g^gy g^t; shown in Fig. 2 of Ref. |^ gives just the ratio between the 
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PSD including the FGM energy dependence of these effective s.p.l. densities, and the ESM state 
densities with only the finite- well depth correction. Both kinds of PLDs are shown in Figs. 9(a) 
(Cases 18/18A) and 9(b) (Cases 18B/18C) for either the normal Fermi energy Fq or the average 
effective value Fi=lA MeV. The analysis concerns the exciton configurations corresponding to 
the first two-body interaction in nucleon-induced PE reactions, which are the most important for 
PE-reaction calculations. 

The two-component PSD formula taking into account the surface effects for the initial target- 



projectile interaction was given as 



PT^.tii^.py.ti^.yn Lj. 



(65) 



where 



A{p^, K,p,, z,^) = ^ + ^ _ Pl + K + n^ _ Pl + hl + n, ^gg^ 

gn gu 4:g^ ^gy 

is the simpler form of the Pauli correction for the Fermi level placed halfway between the last 
occupied s.p.l. in the ground state of the nucleus and the first vacant s.p.l. The same formalism 
as in the one-fermion system case has been used for the average-excitation energies for excited 
particles and holes. The only additional assumption (^6]) concerns the definition of the average 
excitation energies E^^ and E^, for protons and neutrons respectively. 

The PLDs for the two-fermion system exciton configurations (p^, /i^, p,^, /ii,) = (2100), (1110), 
(1100), and (1001) are also compared in Figs. 9(c) (Cases 18A/19/19A) and 9(d) (Cases 
18C/20/20A) with the one-fermion PSDs of the related configurations (p, h) = {21) and (11), as well 
as with the results obtained by using the method involving the two-fermion system correction. A 
behavior specific of the particle-hole bound state density arises for the two-fermion-system exciton 
configurations with only one hole, due to the use of Eq. 

2.5. Composite (recommended) PSD formulas 

The completeness of a composite PSD formula can be obtained by the inclusion and optional 
use of the various ESM corrections : 



extension to the case of the finite well depth and bound states [p^,p!9|,p2 



advanced pairing correction by Fu [T^ and Kalbach |]T6[ added to the Pauli correction of 



Williams, proved to be still in good agreement with exact calculations [|r^-[TP|; 

- inclusion of the shell effects together with the pairing correction [jl5[ and use of the usual 
level-density parameters also for the PSD/PLD calculation; 

- energy dependence of the single-particle state densities as well as inclusion of the surface 
effects in the case of the first two steps of multistep processes ||2^ , which appear to be most 



significant for the exciton configurations mainly involved in the description of PE reactions. 
Therefore, in the one-fermion system case, a composite (recommended) PSD formula is 

p\h\[n — 1)! 

where 
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9p{p,h) = giF + Up) (68a) 

gh{p,h)=g{F-Uh), (68b) 
with the average excitation energies for particles and holes 

nUp,h,E,F) ^''^^ 

u, = , (69b) 



h 



where 



i=0 jr=0 



X e(E - Ethresh -S-iB- jF), (70) 



and 



On the other hand, the two-component partial state density has the general expression 

uip^,k^,p^,k^,E)- pjhjpjhjin-iy. hhhh^ ^ 

xCll Cil C;: Ct [E - Ak{p^, K,p,, K)-S- i^B^ - j^F^ - i,B, - j^FX'' 

X e{E - Etresh -S- i^B^ - j^F^ - kB, - j,F,), (72) 

where 

AxiPTT, K,Pu, K) = AxiPn, K) + Ak{Pu, K), (73) 



the two-fermion systems being considered at the average excitation energies (PB|), respectively. 
The average s.p.l. densities in Eq. (^) are evaluated within the same distinct fermion systems. 
One may observe some difference between the one- and two-component formulas concerning the 
calculation of the nuclear potential finite-depth correction, following the original formalisms. How- 
ever, a two-fermion PSD formula quite similar to the one-component case and consistent with the 
average excitation energies (^6|) may be obtained with minor changes. 
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The optional provision of the advanced corrections to the ESM formulas by Eqs. ( |67D and 
([7^) , respectively, is illustrated in Fig. 10. Thus, no correction to the Williams formulas for either 
the one-fermion system. Fig. 10(a) (Cases 21/21A), or two-fermion system. Fig. 10(b) (Cases 
22/22A), determines PLD values obtained with the composite formulas which are so close or even 
identical to the original ones, for the simplest configurations Iplh. The same is true for the bound- 
state case with only the finite-well depth correction, the comparison with the predictions of the 
Oblozinsky formulas fl^ being shown for specific configurations as well as the total state density 
given by the corresponding PLDs sum. These results are given in Figs. 10(c) (Cases 23/23A) 
and 10(e) (Cases 23A/24A), respectively, for the one-fermion system, and in Figs. 11(a) (Cases 
23B/24B/25B) and 11(c) (Cases 23B/24B) for the two-fermion system formulas. 

The specific behavior following the inclusion of the average energy- dependent formalism is 
pointed out by comparison with the predictions of the Oblozinsky formulas for particular 
configurations and the total state density given by the PLDs sum. This is shown in Figs. 10(d) 
(Cases 23/23C) and 10(f) (Cases 23C/24C), respectively, for the one-fermion system, and Figs. 
11(b) (Cases 23C/24C/25B) and 11(d) (Cases 23D/24D) for the two-fermion system formulas. 
The enhancement of the "bound-state effect" (i.e. reaching a PSD maximum value followed by 
vanishing at a certain higher-energy limit) due to the addition of the energy-dependent s.p.l. 
densities is obvious in Figs. 10(d) and 11(b). A particular feature is shown comparatively in Figs. 
ll(a),(b) for the one- and two-fermion system configurations with Iplh, namely the effects of the 
sequential releasing of the bound-state and finite-well depth conditions. 

The comparison between the predictions of the composite formula, using the FGM energy- 
dependent s.p.l. densities, and the ESM formula including the advanced pairing correction [p!^ , p!^] 
is shown in Fig. 12(a) (Cases 26/26A). The same but for the ESM one- and two-fermion system 



formulas including the exact calculation of the Pauli correction |[T7|-[T9| are shown in Figs. 12(c) 
(Case 28) and 12(d) (Case 28B). The corresponding bound-state PLDs are identical as long as the 
average excitation energies of the two kinds of excitons are lower than the respective limits, while 
next the "bound-state effect" is well increased within the average energy-dependent formalism. 

Finally, the changes due to the inclusion of the surface effects for exciton configurations with 
one and two holes, are illustrated in the case of both fermion systems in Figs. 12(a) and Fig. 
12(b) (Cases 27/27B). These changes are higher than, and increasing, the ones due to the energy- 



dependent s.p.l. densities. The comparison with the original results |2^ not including the pairing 
correction in Fig. 12(b) points out also that taking into account the pairing effects is insignificant 
at higher energies but dominant in the threshold region. 



2. 6. The partial level density 

2.6.1. One- component Fermi gas formula 

The level density of ]9-particle-/i-hole configurations of excitation energy E and nuclear spin 
J, in the one-fermion formulation, is 



p(p, h, E, J) = uj{p, h, E) R{n, E, J) 
where the spin-distribution formula has the general expression i 



R{n,E,J) 



2J 



2{27Ty/^a%E,n 



-exp 



-(^ + 1/2)^ 
2a2(E,n) 



(74) 



(75) 



and the spin-cutoff factor a'^{E,n) may have the simplified energy- and exciton-configuration 
dependence 
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a\E,n) = (ln4) ( — 
Krir 



E-Ui 



th 



E 



where, in addition to the quantities from the previous section, 

X = -0.413 + imin/n^y/^ - 0.226(n/nc 



and is the spin-cutoff factor at the critical thermodynamic temperature Tc 
a1 = gTc{m?). The differences between Eq. (|76D and the spin-cutoff factor a] 



(76) 
(77) 

2A0/3.5, namely 
in) ~ 0.16nA2/3 



initially introduced by Feshbach et al. |^ concern firstly the threshold existing in the former case. 
Secondly, the result of Feshbach et al. is more appropriate for exciton configurations around the 
most-probable exciton number n, which are important to the compound-nucleus contribution; at 
the same time the larger value of about 0.28nA^/3 obtained by Reffo and Hermann [^] is better 
for the n=2 component - the PE-dominant contribution. However, a'^{E,n) can be used for all 
exciton numbers since the respective values are closer to the results of Reffo and Hermann around 
n=2, and also in good agreement with ap{n). 

The sum over the nuclear spins, of the PLDs given by Eq. (0), corresponds to the partial 
level density p{p, h, E) related to the PSD of the same configuration by means of the closed formula 



p{p, h, E) 



{2ny/^a{E,n) 



(7J 



On the other hand, the sum of the PLDs over all allowed particle-hole numbers p = h should 
provide the usual nuclear level density involved in the statistical model calculations of compound- 
nucleus processes, which has also the one-fermion ESM formula 



Pi{E,J)=u;,{E)R{E,J), 



with the spin-distribution 

R{E, J) 



2J+1 



2(27r)i/2^3(^) 



exp 



2a\E) 



(79) 



^0) 



For the spin-cutoff factor cr'^{E) in the above equation, which should be consistent with the value 
<Jh{E) used in compound-nucleus (Hauser- Feshbach) calculations, the average of a'^{E,n) over n 
was adopted 



a\E) 



Ep=h=i ^(P, h, E) 



where the PSD uj{j), h, E) may be given by, e.g., Eqs. ( pT]) or 
to be consistent with the general form of a'jj{E) 



al{E) = g{m^)t 



(81) 

This average was also proved 

(82) 



where the value (m^) =0.24^4^/'^ was assumed [p2 



Furthermore, the sum of the PLDs over both p (with the restriction p=h) and spins gives the 
total level density 



Pi{E)= E T.p{P^h,E,J) 
p=h=i J 



13) 



which is related to the one-fermion ESM total state density by 
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Pi{E) 



{2Txy/^a{E) 



The comparison of the values given by the above closed formula and the sum (^) (i.e. Eqs. (52) 
and (53) respectively of Ref. [|32|) for the excited nucleus ^^Ca is shown in Fig. 13(a) (Case 11). 
The partial and total state and level densities corresponding to the two-component ESM 



formulas can be obtained by using either the renormalization method p2| , |3^ or the adjustment 
method of the two-fermion formula parameters [15,32|. The former, illustrated in Fig. 4(d), offers 
the advantage of using only one set of e.g. BSFG model parameters for both the PLDs involved 
in PE calculations and the nuclear level densities for Hauser-Feshbach calculations. 



2.6.2. Two-component Fermi gas formula 



The two-fermion level density formula is [0 

p(p7r, Pu, K, E, J) = uj{pn, K, Pu, K, E, J) R{n^, n^, E, J) 
where the spin-distribution formula has the similar general expression 



R{n^,n^y,E, J) 



2J+1 



2{27ry /^a^{E,n^,n, 



-exp 



-(^+1/2)^ 

, ?T"7i-, Tip 



i5) 



(86) 



and the spin-cutoff factor (j'^{E, n.^, Uy) for two kinds of fermions is defined as the sum of the two 
one-fermion components 

cr^(i?, n^, Tiy) = cr^{E^, n,r) + (^'^{E^, n^) . (87) 
The mean-gap approximation and the approximations (|4^) for i?^ and Ep as well as the parame- 



terized function for the one-fermion system |^ are used in this respect. 

The sum of the two-fermion PLDs over spins and both and p^ (with the restriction p.^ = 
and Pp = hp) provides the total level density 

P2{E) = Yp{p^,K,Py,K,E,J) (88) 

p^=hTr,Pv=hv J 

which is related to the two-fermion ESM total state density by 

0J2{E) 



P2{E) 



(27r)i/2ff(E) • 



(89) 



The same averages of a'^{E,n) over n as for the one-fermion formulas are used for the above 
spin-cutoff factor 



a\E) 



Ep^=h^,p,=K ^{P^^ Pu, K, E) a'^{E, 

Hp^=hTr,Pu=hi, ^{Pn, hj,, Pu, hp, E) 



n] 



(90) 



where a;(p^, hT,,Pu, hp, E) is given by Eq. 

The comparison of the values given by the above closed formula and the sum (|8^) for the 
excited nucleus ^^Ca is shown in Fig. 13(b) (Case IIB). Fig. 14 (Cases 11/llB) shows the compar- 
ison between the average spin-cutoff values cT|(n=2) = [cr|(i?, 727^=2, n,^=0) + a2{E,nT,=Q,np=2)]/2 



36| and cr^(n = 2) corresponding to the one-fermion system, as well as of the average values (|9 



noted (T2 and the corresponding al (both of them noted "all n" in Figs. 13 and 14). The related 
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PLD and total level density values are shown in Fig. 13(c) (Cases 11/llB) while the comparison 
between the two-fermion system PLD values and the ones obtained by using the two-fermion cor- 



rection |2|,g is given in Fig. 13(d) (Cases IIB/IIE). 



3. Program organization 

The program PLD. FOR is the collection of the above-described algorithms developed until 
now and widely used for PSD/PLD within nuclear model calculations. It also includes the rec- 
ommended (combined) PLD formulas given in this work. Fourteen different PSD formulas are 
available as FORTRAN?? functions to be used as they are in various applications or codes. The 
one method involved in calculating the PLD is given within a subroutine only for the tabular print- 
ing of the results, while the replacement of the SUBROUTINE statement by the FUNCTION one 
is immediate. 

Second, the PLD. FOR has been organized so that various formulas and versions may be 
tried as well as the comparison between their predictions. This could also be useful for further 
development of the PSD calculation methods. 

Since the first aim of this work has been to provide tools for users of PSD/PLD, the opti- 
mization of the respective procedures is made firstly in this respect. The possibility of using these 
functions independently has the related drawback of increasing the execution time. A proper use, 
e.g., of the PLD. FOR for PSD calculation with exact Pauli-correction term of Mao Ming De and 



Guo Hua |19] would involve the calculation of the exact coefficients B{p,h,X) only once in the 



main program. One should keep this aspect in mind if PLD. FOR will be used on low-speed PCs. 
3.1. Subprograms 

The MAIN program reads the input data which are listed and described in Table 1 in reading 
order (including the names of variables which are also used below and stand for various quantities 
given in the previous section). The formatted read is used for all data just to make possible the 
input of only few of them, while by-default values may be involved for the rest. Then, the partial 
state densities w{p,h,E) are calculated by using the specified formula, as well as eventually the 
related partial level densities D{p, h, E, J). In the general case when the calculation is carried on 
for all pairs of the exciton numbers p=h, it is also done for (i) the nuclear state density w{E), 
as the PSD-sum over all allowed exciton numbers for which the PSD is higher than the value 
WMINACC=0.1 MeV-\ or (ii) the total level density D{E) as the PLD-sum over the exciton 
numbers and the nuclear angular momentum J. At the same time the corresponding values 
Wasym{E) or Dasym{E), respectively, are calculated by means of the ESM closed formulas, in 



order to make possible a test of the overall consistency. The two-fermion system correction is 
involved optionally when the one-component Fermi gas formulas are used. 

The subroutine PRINTIN prints the type of the PSD/PLD formula and the parameters 
used in calculation. 

The subroutine PRINTWN tabulates the calculated values w{p,h,E) or D{p,h,E) (the 
latter being the sum over J of D{p, h, E, J), i.e. the total level density for a given exciton config- 
uration), for either (i) some given particle-hole configurations, or (ii) all pairs of equal numbers of 
excited particles and holes. In the latter case the values w{E) or D{E) are also printed within a 
first table including the PSD/PLDs for p=h from 1 to ?, while the corresponding closed- formula 
values Wasym{E) or Dasym{E) are given within the second table including the PSD/PLDs for 
p=h from 8 to 16. A table of the values D{p, h, E, J), but only for the last excitation energy 
E involved in one calculation, is also printed in the case of the PLD calculation for a particular 
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exciton configuration. 

The functions WILl and WIL2 calculate the partial state density w{p, h, E) for a given exci- 
ton configuration, by using the Williams [0] one- and two-fermion system formulas, respectively. 

The functions WOBl and WOB2 calculate the partial state density for a given particle-hole 
configuration by means of the Betak-Dobes one- and two-fermion system formulas, respectively 



12| , with the nuclear potential finite-depth correction. The calculation of the bound-state density 
according to the Oblozinsky formulas |]14| is carried out if a value is specified for the nucleon 
binding energy. 

The functions WFUl and WFU2 calculate the PSD for a given exciton configuration by 
means of the one- and two-fermion system formulas, respectively, including the advanced pairing 



correction by Fu |]T5|. 



The functions WKl and WK2 calculate the PSD for a given exciton configuration by using 



the improved implementation of the pairing correction by Kalbach [|l^] within the one- and two- 
fermion system formulas, respectively. 

The functions WK3 and WK4 calculate the PSD for a given exciton configuration by using 
the FGM energy-dependence of the single-excited particle and single-hole state densities, and/or 



the finite-depth correction including the nuclear-surface effects introduced by Kalbach ||2^ within 
the one- and two-fermion system formulas, respectively. 

The functions WMl and WM2 calculate the PSD for a given exciton configuration by using 
the exact calculation of the Pauli-exclusion effect and the pairing correction by Kalbach [|l| 



within the one- and two-fermion system formulas, respectively. 

The functions WRl and WR2 calculate the PSD for a given exciton configuration by using 
the composite (recommended) formalism including optionally (i) the improved implementation of 
the pairing correction by Kalbach , (ii) the FGM energy-dependence of the single-excited par- 
ticle and single-hole state densities, and/or (iii) the finite-depth correction including the nuclear- 
surface effects introduced by Kalbach ^2[, within the one- and two-fermion system formulas, 
respectively. 

The function PFU calculates the advanced pairing correction by Fu [|1^, for the one- 
component Fermi-gas. 



The function AK calculates the advanced pairing correction by Kalbach |]T6[, for the one- 
component Fermi-gas. 

The function FDCO calculates the nuclear potential finite-depth correction factor 
f{p, h, E, F) p2[ for the one-component Fermi-gas. 

The function FDC calculates the nuclear potential finite-depth correction factor f{p + 
1, h, E, F) but for the case of bound states. 

The subroutine SUBPLD calculates the partial level density D{p, h, E, J) for a given particle- 
hole configuration, excitation energy E and nuclear spin J, as well as the respective total level 
density D{p,h,E) as their sum over J, by using the partial state density w{p,h,E) and the 
formalism of Fu . 

The function SIG2FU calculates the spin-cutoff factor for a given excited particle-hole con- 
figuration . 



The function FCTR calculates the factorial of natural numbers. 

4. An illustrative test run 

The sample case 23C (see also Fig. 10) is here discussed since it documents few of the specific 
features of the program PLD at once. Thus, the particle-hole bound state densities for some of the 
few-exciton configurations analyzed by Oblozinsky (Fig. 1 of WM) are calculated by using both 
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the composite formula and the one of Oblozinsky. In the former case the FGM energy dependence 
for the s.p.l. density is additionally taken into account, while the input data correspond to the 
reference work The comparison of the results obtained with the two formulas is possible 

within the same table with the columns of results given in the input-data reading order (the table 
in the attached output copy below is reduced to the first half of the involved excitation energies). 
The corresponding curves are shown in Fig. 10(d). One note should concern the printed type of 
the PSD formula used: it denotes only the last one when the option parameter IC0NT=2 is used 
and more than one formula are involved in calculation. 

An additional calculation is included in this case with respect to its title. To make possible 
the comparison between the PSD for a given exciton configuration with equal numbers of holes 
and excited particles which is calculated both specifically and within the general case for all pairs, 
the latter calculation is also made. Moreover, in the attached copy of the reduced output only the 
first two tables for this calculation are shown (which would correspond to the use of the option 
parameter value IC0NT=-1). One may thus find the printed values of the total state density w{E) 
obtained as the PSD-sum, in the first table, and the related Wasym{E) closed-formula values, in 
the second table. The corresponding curves are shown in Fig. 10(e). However, this part of the 
case output consists of four tables which correspond to the maximum numbers p=h=25 for which 
the PSD value at higher limit of excitation energies is yet higher than WMINACC. 

Actually, the PSD/PLD values for exciton numbers higher than n are usually of less interest, 
while the PLD-program output includes them just for the sake of completeness. In the possible 
case that an output table would include only zero values, it is omitted from print. We may add 
that questions may arise for calculations for all pairs p=h when quite large excitation energies and 
the two-fermion system formulas are involved, due to the limits of vectors. However, the use of 
the respective FUNCTION subprograms for given exciton numbers and energy, i.e. the usual case 
of nuclear reaction cross-section calculations, is always straightforward. 
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TEST RUN (REDUCED) OUTPUT 



PARTIAL STATE/LEVEL DENSITIES CALCULATED FROM THE FOLLOWING PARAMETERS 

Case 23C: C-f ormula(g-FGM)/0blozinsky(1986) : one-f . bound PSD for configs. 

NUCLEUS: CHARGE NO. Z= 0. MASS NO. A= 0. PAIRING Up= .000 MEV 

ONE-FERMION FORMULA: 

P.OBLOZINSKY, NUCL.PHYS. A453, 127(1986) , Eqs.(7,9) 

SINGLE-PARTICLE STATE DENSITY: G= 8.000 /MEV 

FERMI ENERGY: F= 32.000 MEV 

NUCLEON BINDING ENERGY: B= 8.000 MEV 
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PARTIAL STATE/LEVEL DENSITIES CALCULATED FROM THE FOLLOWING PARAMETERS 



Case 23C: C-f ormula(g-FGM)/0blozinsky(1986) : one-f . bound PSD for configs. 

NUCLEUS: CHARGE NO. Z= 0. MASS NO. A= 0. PAIRING Up= .000 MEV 

ONE-FERMION FORMULA: 
COMPOSITE FORMULA 

SINGLE-PARTICLE STATE DENSITY: G= 8.000 /MEV 
FERMI ENERGY: F= 32.000 MEV 

AV. EFF. FERMI ENERGY/ 1ST-2ND STEPS: Fl= 32.000MEV 
NUCLEON BINDING ENERGY: B= 8.000 MEV 
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15, 


.00 


.154E+11 


. 660E+09 


.313E+08 


. 194E+06 


.000 


.000 


.000 


.000 
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16, 


.00 


.362E+11 


.247E+10 


. 176E+09 


. 225E+07 


.112E+04 


.000 


.000 


.000 


17, 


.00 


.833E+11 


.826E+10 


. 827E+09 


. 185E+08 


. 348E+05 


.000 


.000 


.000 


18, 


.00 


. 187E+12 


.252E+11 


.338E+10 


.118E+09 


. 554E+06 


.000 


.000 


.000 


19, 


,00 


.412EH-12 


.707E+11 


.123E+11 


. 628E+09 


. 577E+07 


.191E+04 


.000 


.000 


20. 


.00 


.890E+12 


. 186E+12 


.407E+11 


.286E+10 


.446E+08 


. 538E+05 


.000 


.000 


21. 


.00 


. 189E+13 


.458E+12 


. 124E+12 


.115E+11 


. 276E+09 


.822E+06 


.000 


.000 


22, 


.00 


.394E+13 


.107E+13 


.351E+12 


.417E+11 


. 144E+10 


. 843E+07 


. 126E+04 


.000 


23, 


.00 


.811E+13 


.240E+13 


.933E+12 


.138E+12 


.654E+10 


.652E+08 


.390E+05 


.000 


24, 


.00 


. 164E+14 


.514E+13 


.235E+13 


.426E+12 


.265E+11 


.408E+09 


. 644E+06 


.000 


25. 


.00 


.328E+14 


. 106E+14 


.562E+13 


. 122E+13 


.971E+11 


.216E+10 


.705E+07 


287. 


26. 


,00 


. 647E+14 


.210E+14 


. 129E+14 


.332E+13 


.328E+12 


.lOOE+11 


. 577E+08 


. 124E+05 


27, 


.00 


.126E+15 


.405E+14 


.283E+14 


.855E+13 


. 103E+13 


.416E+11 


. 380E+09 


.253E+06 


28, 


.00 


.243E+15 


.757E+14 


.601E+14 


.210E+14 


.303E+13 


.157E+12 


.211E+10 


.324E+07 


29. 


.00 


.462E+15 


.138E+15 


. 123E+15 


.494E+14 


.845E+13 


.545E+12 


.102E+11 


.299E+08 


30. 


,00 


.872E+15 


. 244E+15 


.245E+15 


.112E+15 


.224E+14 


. 176E+13 


.443E+11 


.217E+09 


31. 


,00 


. 163E+16 


.423E+15 


.473E+15 


. 244E+15 


.565E+14 


. 536E+13 


. 174E+12 


.131E+10 


32, 


.00 


.301E+16 


.716E+15 


.889EH-15 


.517E+15 


. 137E+15 


. 154E+14 


.631E+12 


.683E+10 


33, 


.00 


.551E+16 


.119E+16 


. 163E+16 


.106E+16 


.320E+15 


.421E+14 


.213E+13 


.315E+11 


34, 


.00 


. lOOE+17 


.193E+16 


.292E+16 


.211E+16 


.721E+15 


.llOE+15 


.672E+13 


.132E+12 


35. 


,00 


. 180E+17 


. 308E+16 


.512E+16 


.411E+16 


. 157E+16 


.276E+15 


.201E+14 


. 504E+12 


36. 


,00 


.322E+17 


.483E+16 


.879E+16 


.778E+16 


.333E+16 


. 665E+15 


.571E+14 


. 179E+13 


37. 


.00 


.572E+17 


.745E+16 


. 148E+17 


. 144E+17 


.687E+16 


.155E+16 


. 155E+15 


.594E+13 


38, 


.00 


.lOlE+18 


.113E+17 


.245E+17 


.261E+17 


. 138E+17 


.350E+16 


.403E+15 


.186E+14 


39, 


.00 


. 176E+18 


.169E+17 


.398E+17 


.465E+17 


.270E+17 


.767E+16 


.lOlE+16 


.553E+14 


40. 


.00 


.306E+18 


. 250E+17 


. 637E+17 


.809E+17 


.517E+17 


. 163E+17 


. 244E+16 


. 157E+15 
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Figure Captions 



FIG. 1. Particle-hole state densities for the given ph configurations, and the sum of the PSD for all 
allowed exciton numbers p=h (crosses), obtained with (a) the Williams' one- and (b)-(d) 
two-fermion formulas for (a) the generic value g=l MeV^\ (b) g^=g^=g/2=l MeV~\ and 

(c) ,(d) the nucleus ^^Nb and the phenomenological value g=A/13 MeV^^. There are also 
shown the total nuclear state-densities given by the ESM formulas for (a) the one- and (b)- 

(d) two-fermion systems [0 (solid curves), compared with (c),(d) the one-fermion system 



closed formula (dotted curves) and (d) the sum of the renormalized one-fermion PSDs [0 
(dashed curve). 

FIG. 2. Particle-hole state densities for the 2plh and 3p2h configurations, obtained with the ESM 
one-fermion formula with (dashed and dotted curves, for various F values) and without 
(solid curve) the potential finite-depth correction |T2|, for the value g=l MeV~^. 



FIG. 3. Particle-hole bound state densities obtained with the Williams' one-fermion formula includ- 
ing the nuclear-potential finite depth [|I^] for (a) the given ph configurations (solid curves) 



and (a),(b) the configuration Iplh, releasing consecutively the finite-depth potential (dashed 
curve) and the bound state conditions (dotted curve) by means of large F- and S-values, 
(b),(c) their sum for all allowed exciton numbers p=h without/with the two corrections 
(upper/lower crosses), and (c) the similar sum of the two-fermion PSDs (upper/lower x) 
compared with the one- and two-fermion system closed formulas (dashed and solid curves, 
respectively), for the values g=8 MeV'^ F=32 MeV, and B= 8 MeV. 
FIG. 4. PSDs obtained with the Williams one-fermion formula (a) without/with the advanced pairing 
correction [^] (dashed/solid curves) and (b) above the threshold energy for each configu- 
ration (dotted curves), leading to the sum of the PSD for all allowed exciton numbers p=h 
(solid curve), compared with the closed formula (dashed curve) for the values g=14: MeV~^ 
and Ao=l MeV, as well as (c) for g=4: MeV~^ and constant pairing correction Up of 0, 2, and 
4 MeV, and (d) for the one/two-fermion system formulas (upper /lower respective curves) 
and the BSFG parameters a=4.12 MeV~^ and A=-1.07 MeV for the excited nucleus ^^Ca 



FIG. 5. Comparison of the PSDs obtained with the Williams one-fermion formula, for the given 
exciton configurations, and the advanced pairing correction of (a) Fu |T5| (dashed curves) 



and Kalbach above the threshold energy for each configuration (solid curves), as well 
as of (c) the latter and the PSDs with exact calculation of the Pauli-principle correction 
|[T9[| (dotted curves), (b) the sum of the Kalbach PSDs (shown also below the threshold) for 
all allowed exciton numbers p=h (solid curve) compared with the closed formula (dashed 
curve), for the values (7=14 MeV~^ and Ao=l MeV, and (d) the same but for the exact 



formula and the value g=8 MeV ||19|| . 



FIG. 6. Comparison of the particle-hole bound state densities for the given ph configurations and 
the values g=8 MeV~\ F=32 MeV, and B=8 MeV, obtained with the Williams one-fermion 
formula (a) including the nuclear-potential finite depth (dotted curves) and the exact 
calculation of the Pauli-principle correction |]19| (solid curves), as well as with (b) the latter 
form but without /with the advanced pairing correction and parameter value Ao=l MeV 
(dotted/sohd curves). 

FIG. 7. The sum (solid curves) of (a) the two-fermion system PSDs ||16|| (dotted curves) for al- 
lowed exciton numbers Pn=hTt and Pu=hu, compared with (b) the closed formula of the 
two-component ESM total state density for the excitation energy decreased by the effective 
pairing correction given by the term P!^_^^^^^f (dotted curve), and (b),(d) the term P7r+j/,e// 
(dashed curve), (c) the constant Up (long-dashed curve), the sum Pf,ff{E, C.„)+xPeff{E, C^) 
with x=l (dashed curve) as well as between and 1 for E/C varying from 1 to 2 (dotted 
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curve), and (d) the Kalbach term Pf,ff{E, C), also compared with the sum of the renormal- 



ized one-fermion PSDs p4[ (including the advanced pairing correction for all allowed 
exciton numbers p=h (long-dashed curve), for the values (7=14 MeV~^ and Ao=l MeV. 
FIG. 8. The sum (solid curves) (a) of the two-fermion system PSDs (dotted curves) for allowed 
exciton numbers p^=/;,^ and Pu=hy, compared with the closed formula of the two-component 
ESM total state density for the excitation energy decreased by the effective pairing correc- 
tion given by (b),(c) the Kalbach effective-energy shift Peff{E,C) (dotted curve), the term 
PiT+u,eff (dashed curve), and (c) the constant shift Up. (d) The comparison of the PSD 



values of Refs. [|T9l (dashed curves) and |T6|, for the given exciton configurations. The global 
parameter values used are Ao=l MeV, and (a),(b) g=S MeV~^ |T^JT9[] and (c),(d) 5'=14 
MeV^i EH. 

FIG. 9. PSDs for (a),(b) the given ph configurations, obtained with the ESM one-fermion system 
formula with (solid dotted) and without (dotted curves) inclusion of the FGM energy de- 
pendence of the s.p.l. density, and (c),(d) the given p-KhT,pyhu configurations, obtained with 
the two-fermion system formula including the FGM energy-dependent s.p.l. density (solid 
curves) compared with the related results of the one-fermion system formula (dashed curves) 
and including the two-fermion system correction (dotted curves), as well as (a),(c) without 
and (b),(d) with the surface effect taken into account by using the average effective well 



depth Fi=14 MeV for the s.p.l. density value 5'o=14 MeV and Fermi energy Fo=38 
MeV. 

FIG. 10. Comparison of the PSD values for the given ph configurations and the sum providing the to- 
tal state density, obtained by means of the given parameter values and (a),(b) the Williams 
(dotted curves) and the composite formula with no additional corrections (solid curves), 
for the one- and two-fermion systems respectively, and (c),(e) the Oblozinsky one-fermion 
formula [Q (dotted curves) and the composite formula including the finite-well depth and 
bound-state conditions (solid curves) as well as (d),(f) the average energy-dependent formal- 
ism. 

FIG. 11. Comparison of the PSD values for the given ph configurations and the sum providing the to- 
tal state density, obtained by means of the given parameter values and (a),(c) the Oblozinsky 
one- and two-fermion system formulas (dotted curves) and the composite formula in- 
cluding the finite-well depth and bound-state conditions (solid curves) as well as additionally 
(b),(d) the energy-dependent s.p.l. densities within the composite formula. The PSDs for 
configurations including Iplh which are obtained as noted in (a),(b) by releasing consec- 
utively the finite-depth potential and bound-state conditions by means of large F- and 
5- values are also shown. 

FIG. 12. Comparison of the PSD values for the given exciton configurations, obtained by means of 
the given parameter values and the composite formula including the finite-well depth, the 
energy-dependent s.p.l. densities (solid curves) and the surface effects taken into account 
by using the effective well depth Fi=14 MeV (dashed curves), and (a) the Kalbach one- 
fermion system formula [|l^] (dotted curves), (b) the Kalbach two-fermion system formula 
(dotted curves), as well as (c),(d) for the bound states and the exact calculation of the 
Pauli-principle correction [1^ (dotted curves). 

FIG. 13. Comparison of the total level density given by the sum of the ESM (a) one-fermion and (b) 
two-fermion system PLDs over both excited-particle number (s), equal to hole number (s), 
and spins (solid curves) and the related closed-formulas (dashed curves), as well as of the 
closed-formula values obtained by using (c) the two-fermion system average spin-cutoff values 
(j|(n = 2) and the correspondent <y'l{n = 2) for the one-fermion PLD (dotted curves) 
and the average values cr| and af for the related total level densities (solid curves noted "all 
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n"), and (d) the two-fermion system formula (solid curve) and the one-fermion formula with 
the two-fermion correction |]32| , |3^ (long-dashed curve), for the excited nucleus ^^Ca and the 
given parameter values. The corresponding PLD values are shown for configurations with 
two excitons (dotted curves). 
FIG. 14. Comparison of the two-fermion system average spin-cutoff values (T2{n = 2) and the 
corresponding af{n = 2) for the one-fermion PLD (dotted curves), and the average values 
(t| and af (solid curves noted "all n"), for the excited nucleus ^^Ca and the given parameter 
values. 
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Table 1. Sample cases for partial state density (PSD) and partial level density (PLD) calculation 

up to excitation energy E^ax- 



Case 


TITLE (variable on first input record) 




Execution 


No. 




(MeV) 


time" (s) 


1 


Williams(1971)/i' ig.2(one-iermion, asymptotical lorm m 2nd table!) 


60 


0.2 


lA 


Williams (1971) one-termion tormula, tor distinct conngurations 


40 


0.1 


2 


Williams(1971),as rig.2 but two-tcrmion tormula/asymptotical lorm 


40 


0.2 


2A 


Williams(1971) two-lermion formula, lor distinct conngurations 


40 


0.1 


3 


Williams(1971) one-fermion formula/ asympt. form for JNb, g=A/13 


35 


0.1 


3A 


TTT'11* f -\ f \ \ j P * P 1/ I C P"\T1 A/"10 

Williams(1971) two-lermion tormula/ asympt. lorm lor iNb, g=A/13 


35 


0.3 


3B 


Williams(1971) one-fermion+TFC'' formula/asympt.form for Nb, g=A/13 


35 


0.1 


4 


Betak+(197o)/l' ig.l:one-lermion state density + nnite depth corr. 


400 


0.3 


5 


(Jblozinsky(198o)/l'ig.l:one-lermion bound state dens, lor conngs. 


80 


0.1 


5A 


Oblozinsky(1986)/Fig.l: one-fermion state density + finite depth 


80 


0.1 


5B 


Oblozinsky(1986)/Fig.l: one-fermion state density formula 


80 


0.1 


6 


(jblozinsky(198o) onc-lermion bound-state densities: g=8 


80 


0.4 


6A 


Oblozinsky(1986) one-fermion state density + finite depth corr. 


80 


0.4 


6B 


Oblozinsky(1986) one-fermion formula/asymptotical form: g=8 


80 


0.3 


7 


Oblozinsky(1986) two-fermion bound-state densities: g=8 


80 


9.2 


7A 


11 • 1 / \j {* • jj 1 •! |P*"J 1 jl C\ 

Oblozinsky(1986) two-fermion state density -|-nnite depth corr., g=8 


80 


7.0 


7B 


(Jblozinsky(198o) two-termion tormula/asymptotical torm: g=8 


80 


6.8 


8 


Fu(1984)/Fig.3: one-fermion state density with/without pairing 


15 


0.3 


9 


Fu(1984)/Fig.6: one-f. closed-formula values for various pairings 


20 


0.3 


10 


1 1 / ^ * y-^ J \ 1 1 1 • /, fl If* 1 "jl ■ -| X 1 ■ 1 J 

Fu(1984)/Fig.7: one/two-f. closed-formula with various BSFG data 


40 


0.2 


lOA 


Pu(1984)/Fig.7: two-f. closed-formula with two-f. correction 


40 


0.1 


11 


1 — 1 ( -y f\ o \ /i — 1 P • j'llll /I IP 1 

lu(198o)/l ig.o:one-termion partial level dens, sum/closed-tormula 


11 


0.1 


UA 


lu(198o) one-termion partial level density, total rLD tor conng. 


11 


0.1 


IIB 


Fu(1986)/Fig.l: two-fermion PLD for 41Ca: sum/closed-formula 


25 


0.2 


lie 


i*u(1989)/i'ig.l: two-fermion PLD for 41Ca connguration 


25 


0.1 


IID 


T — ^ / /~\ /— t /— \ \ / "1 — T • J f ■ "I — \ -f- "I — V (* At / — ^ / / 1 IP 1 

Fu(1989)/Fig.l: two-fermion PLD for 41Ca /sum/ closed-formula 


25 


0.1 


HE 


Fu(1989)/Fig.l: one-fermion PLD -|- TEC: sum/closed-formula 


25 


0.1 


12 


Kalbach(1987)/Eig.3: one-fermion partial state dens, for configs. 


15 


0.1 


13 


Kalbach(1987) one-fermion partial state density/sum/closed form 


15 


0.1 


13A 


Kalbach(1987) two-fermion PbD sum/closed formula,tor g=8/Up=3.5 


15 


0.3 


13B 


Kalbach(1987) one-fermion PSD + TEC /sum/closed form: g=8/Up=3.5 


15 


0.1 


14 


Kalbach(1989)/ lab.l: one-termion partial state dens, tor conngs. 


32 


0.1 


14A 


Kalbach(1989)/ lab.lb:onc-termion partial state dens, tor conngs. 


128 


0.1 


14B 


Kalbach(1989)/Tab.lc:one-fermion partial state dens, for configs. 


256 


0.1 


15 


Mao Mmg De(1993)/l*ig.l:one-fermion bound PSD -|- exact Pauli corr. 


100 


106.7 


15A 


Mao Mmg De(1993)/Eig.2/Oblozinsky(86):bound PSD-Fexact Pauh corr. 


100 


42.9 


16 


Mao Ming De(1993)/Eig.3:bound PSD+exact Pauli corr.+pairing corr. 


100 


64.1 


17 


Mao Ming De(1993):onc-fermion PSDs+pairing corr. /sum/closed form 


15 


123.4 


17A 


Mao Ming De(1993)/Kalbach(87):one-f. PSDs: exact Pauh corr. effect 


15 


134.6 


17B 


Mao Ming De(1993):two-fermion PSDs-|-pairing corr. /sum/ closed form 


15 


351.5 


17C 


Mao Ming De(1993):l-f.PSDs+pair.corr.(2-3.5)+TEC/sum/closed form 


15 


257.4 


17D 


Mao Ming De(1993)/Kalbach(87):two-f. PSDs: exact Pauh corr. effect 


15 


427.1 
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Table 1. - continued 



Case 


iiiijii/ i^vaiiaDie on nibt input lecoiQj 


XT 

^max 


iL/xecuiion 


l\Tr> 
iNO. 




^^lVleV ) 


time" (s) 


1 S 

io 


i\.aiDacn(^iyooj/ (^r ig.zj.one-ierniion rou witn rinite ueptn L/0r..ri — oo 


1 nn 
iUU 


n 1 
U.i 


ISA 
ioA 


ivaiDacni^iyooj/ (^r ig.zj.onc-i. iivnergy-ciep. s.p.s. aensity+r .r — oo 


1 nn 
iUU 


n 1 
U. 1 


ioiJ 


ivaiDacni^iyooj/ (^r ig.zj.one-i. rou witn ruu+suriace eiiects. ri — 14 


1 nn 
iUU 


n 1 
U. 1 




ivaiuacn^^iyoo ) / \" ig-z j.ii/-Qep. b.p.b.uenbity-t-r uv^.r — o(5-rbuii.ei..r i — 14 


1 nn 


n 1 

U. 1 


1 Q 


ivaiDacin^iyoOj/ i^r ig.zj. two-ierniion rou -r ru*^ yr — ooj, g-rvjivi lor 


1 nn 


n 
u.z 


1 Q A 


ivaiDacin^iyooj. one-iennion rou -|- ru»^ — oo) -\- Lrv^, lor conngb. 


1 nn 


n 1 
u.i 


on 


ivaiDacni^iyooj/ (^r ig.zj. two-iermion rLiL) -\- ruu [r i — i^J; g-rvjivi, lor 


1 nn 


n o 
u.z 


on A 


ivaiDacni^iyooj. one-iermion rLiU + ruu yr — i4j + iru 


1 nn 

iUU 


n 1 

U. 1 


Ol 


C-forniula/Willianis(1971)/Fig.2: one-fermion PSD /sum/closed form 


fin 
oU 


n Q 
u.o 


01 A 


C-forniula/Willianis(1971)/Fig.2: one-fermion PSD for configs. 


fin 

DU 


n 1 
u.i 


OO 


C-formula/Williams(1971): two-fermion PSD /sum/closed formula 


/in 

4U 


n 

u.o 


00 A 


C-formula/Williams(1971): two-fermion PSD for distinct configs. 


/in 

4U 


n 1 

U. 1 


OQ 
ZO 


C-formula/Oblozinsky(1986)/Fig.l: one-f. bound PSD for configs. 


en 


n 1 
U. 1 


0*? A 
ZoA 


C-formula/Oblozinsky(1986) ! one-fermion bound PSD/sum/closed form 


sn 


nil 
U. 11 


OQR 


<^-iormuia/ wDlozmsKy^^iyoDj . two-iermion Douna rou/sum/cioseci lorm 


sn 


O/l s 
Z4.0 


ZOv> 


v>-iormuia^^g-r vjivij/ WDiozmbKy^^iyoDj . one-i. DOuna triju lor conngb. 


sn 


o.o 


ZoU 


o-iormuia(^g-r (orivij/ ijDiozmsKyi^iyoOj . z-i. douuq r ou/sum/cioseQ lorm 


sn 

oU 


OS 1 
Zo. 1 


O/l 
Z4 


C-formula/Oblozinsky(1986) : one-fermion PSD -|- FDC for configs. 


sn 
oU 


n 1 
U. 1 


0/1 A 
Z4A 


C-formula/Oblozinsky(1986) : one-f. PSD -|- FDC /sum/closed form 


sn 


1 n 
l.U 


O/l R 


U-iormula/UDlozmsKy(^iyoDj: two-i. rou + rUU /sum/closed. lorm 


sn 


1 7 Q 


o/ir" 

Z'i*^ 


<^-iormuia\^g-r Vjivij/ vjDiozmsKy ^^lysDj . one-i. r ou-|-r uv_>/sum/cioseci lorm 


sn 


i.O 


Z4U 


^^-lormuia^^g-r vjrivij/ UDiozmsKy^^iyooj. xwo-i. r ou-|-riyo/sum/cioseQ lorm 


sn 


10 7 

ly. 1 


oc; 

ZO 


C-formula/Oblozinsky(1986) : one-fermion ESM PSD for configuration 


sn 


n 1 

U. 1 


01^ A 
ZOA 


C-formula/Oblozinsky(1986) ! one-fermion ESM PSD/sum/closed form 


sn 


1 n 

l.U 


ZOJJ 


v^-iormuia/ wDlozmbKy^^iyoDJ . two-iermion Hioivi troiJ /bum/ciobea lorm 


sn 


1 fi 7 
ID. / 


26 


C-form(g-FGM) /Kalbacli(1987):one-f.PSD-FFDC(F=38MeV) /sum/closed form 


80 


2.3 


26A 


C-form(g-FGM) /Kalbach(1987):l-f.PSD+FDC(Fl=14MeV) /sum/closed form 


80 


2.5 


27 


C-form(g-FGM)/Kalbach(1985):two-f.PSD+FDC(F=38MeV)/sum/closed form 


80 


46.7 


27A 


C-form(g-FGM)/Kalbach(1985):2-f.PSD+FDC(Fl=14MeV) /sum/closed form 


80 


46.9 


28 


C-form(g-FGM)/Mao Ming(1993): l-f.PSD+FDC(F=32MeV) /sum/closed form 


80 


276.6 


28A 


C-form(g-FGM) /Mao Ming(1993) :2-f.PSD+FDC(F=32MeV)/sum/elose(l form 


80 


1679.3 



" For calculation carried out on PC Pentium/ 166MHZ Intel. 
Two-fermion system correction (TFC). 
Nuclear potential finite-depth correction (FDC). 
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Table 2. Input data in reading order. 



Record Variable FORMAT Meaning 

1 NE 2I3,74A1 if NE<0, |NE| is the number of excitation energies to be read (<200); 

if NE>0, the excitation energies are (FL0AT(I),I=1,NE) (MeV) 

lOPTJ spin distribution (PLD calculation) option: 

= 0, Lo{p,h,E) or uj{pt^, h-j^jPu, hp, E) are calculated 

= 1, p{p, h, E, J) or p^Ptt, fiT^^Pu, hy,E, J) are calculated 

TITLE title of the problem 

2 (E(I),I=1,-NE) 8F10.5 [record to be read only if NE<0] 

excitation energies (MeV) for PSD/PLDs calculation 

3 IMOD 2I3,7F10.5 option for PSD formula (odd/even: one/two- fermion system formulas): 

=-1, 0: composite (recommended) formula 
= 1, 2: F.C. Williams, Nucl. Phys. A166, 231 (1971) 
= 3, 4: P.Oblozinsky, Nucl. Phys. A453, 127 (1986), Eqs.(7,9) 
= 5, 6: C.Y. Fu, Nucl. Sci. Eng. 86, 344 (1984) 

= 7, 8: C. Kalbach, Nucl.Sci.Eng.95,70(1987), Z.Phys.A 332,157(1989) 
= 9,10: C. Kalbach, Phys. Rev. C 32, 1157 (1985) 
=11,12: Mao Ming De Gua Hua, J. Phys. G 19, 421 (1993) 
option for the two-fermion system correction: 
= 0: no two-fermion system correction 

= 1: J.M.Akkermans, H.Gruppelaar, Z. Phys. A 321,605(1985), Eq. (9) 
excited-nucleus mass number (it may be omitted if GIN>0.) 
excited-nucleus atomic number 

pairing correction based on the odd-even mass differences 
(if UP=-1. and A>0., Eq.(9) of |3| is adopted) 

4 NPO 2I3,7F10.5 [record to be read only for odd IMOD] number p of excited particles 
NHO number h of holes 

(if NPO=NHO=0, calculation is done for all pairs p=h=l, 2,... 
for which PSD/PLD-value is >0.1 MeV'^) 
GIN single-particle state density G (MeV~^): 

if GIN<0., then G=(6/3.1422)*DR(1) is adopted; 
if GIN=0 and A=0, then G=1.0 is adopted; 
if GIN=0 and A>0, then G=A/13.0 is adopted 
FIN Fermi energy F (MeV); if FIN=0, then F=10^ is adopted 

BIN nucleon binding energy B (MeV); if BIN=0, then B=10^ is adopted 

FUN average effective Fermi energy Fl (MeV); if FIN=0, then Fl=10*3; 

if FIN<0, then constant G is used within the WRl and WR2 functions 



ITFC 



A 
Z 

UP 



33 



Table 2. - continued 



Record Variable FORMAT Meaning 



NPO 
NHO 
NPNO 
NHNO 



GIN 

GN 



FIN 

FN 

BIN 
BN 

FUN 



FIN 



4I3,6F10.5 [record to be read only for even IMOD] number of proton excited particles 
number h^^ of proton holes 
number of neutron excited particles 
number /ijy of neutron holes 

(if NPO=NHO=NPNO=NHNO=0, all configurations {p^=h.,r ,Pu=K) 
with N=2, 4, .. and PSD/PLD sum >0.1 MeV'^ are considered) 
single-proton state density G (MeV""'^): 
single-neutron state density GN (MeV~^): 
if GIN<0., then G=Z/A*(6/3.1422)*DR(1) (see record 6) 
and GN=(A-Z)/A*(6/3.1422)*DR(1) are adopted; 
if GIN=0 and A=0, then G=GN=1.0 are adopted; 
if GIN=0 and A>0, then G=Z/13.0 and GN=(A-Z)/13.0 are adopted 
proton Fermi energy F (MeV); 
neutron Fermi energy FN (MeV); 
if FIN=0, then F=FN=10'' are adopted 

proton binding energy B (MeV); if BIN=0, then B=10^ is adopted 

neutron binding energy B (MeV); 

if BIN=0, then B=BN=106 are adopted 

average effective proton Fermi energy Fl (MeV); 

if F1IN=0, it is adopted the value Fl=10'^; 

if F1IN<0, constant G is used within WRl and WR2 functions 

average effective neutron Fermi energy FIN (MeV); 

if F1N=0, it is adopted the value F1N=106; 

if F1N<0, constant GN is used within WRl and WR2 functions 



6 DR(I) 8F10.5 



7 ICONT 213 



lEND 



[record to be read only if GIN<0] BSFG level density parameters |35|: 
DR(1) - nuclear level density parameter a; 

DR(2) - ratio of effective nuclear moment of inertia to rigid-body value; 
DR(3) - shift of the fictive nuclear ground state 
output and recycle option: 

=-1, print first 2 tables of calculated PSD/PLDs (see PRINTWN description); 

resumption according to lEND; 
= 0, print calculated PSD/PLDs; resumption according to lEND; 
= 1, new case by input-data record 1; results printed at once with previous case; 
= 2, new case by input-data record 3 (same energy grid); 
= 3, calculation for other exciton configuration by input-data record 4 or 5; 
= 4, calculation for other BSFG parameter set by input-data record 6 
recycle option: 
= 0, end; 

= 1, new complete case by input-data record 1; 

= 2, new case by input-data record 3; 

= 3, new case for exciton configuration by record 4 or 5; 

= 4, new case for BSFG parameters by record 6 
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